SimulatingMagnetohydrodynamic Instabilitieswith Conservative Perturbed MHDModel Using Discontinuous Galerkin Method

In magnetically confined plasma research, the understandings of small and large perturbations at equilibrium are both critical for plasma controlling and steady state operation. Numerical simulations using original MHD model can hardly give clear picture for small perturbations, while non-conservative perturbed MHD model may break conservation law, and give unphysical results when perturbations grow large after long-time computation. In this paper, we present a nonlinear conservative perturbed MHDmodel by splitting primary variables in original MHD equations into equilibrium part and perturbed part, and apply an approach in the framework of discontinuous Galerkin (DG) spatial discretization for numerical solutions. This enables high resolution of very small perturbations, and also gives satisfactory non-smooth solutions for large perturbations, which are both broadly concerned in magnetically confined plasma research. Numerical examples demonstrate satisfactory performance of the proposed model clearly. For small perturbations, the results have higher resolution comparing with the original MHDmodel; for large perturbations, the non-smooth solutions match well with existing references, confirming reliability of the model for instability investigations in magnetically confined plasma numerical research. AMS subject classifications: 65M60, 65Z05, 85-08